Selective interaction and its effect on collective motion

Plenty of empirical evidence on biological swarms reveal that interaction between individuals is selective. Each individual’s neighbor is selected based on one or more featured factors. Based on the self-propelled model, we develop a general probability neighbor selection framework to study the effect of four typical featured factors (i.e., distance, bearing, orientation change and bearing change). In this work, two common cases are involved to comprehensively analyze the impact of the four featured factors on the collective motion. One is the flocking, the other is the responsivity to stimulus. The impact of different selection strengths of the featured factors on both cases are investigated. The effect of noise on flocking and different stimulus intensities on responsivity to stimulus are analyzed. This study allows us to get the insight of selective interaction and suggests the potential solution to overcome the trade-off between flocking and responsivity quality.


Selective interaction and its effect on collective motion
Zhicheng Zheng 1 , Xiaokang Lei 2 & Xingguang Peng 1* Plenty of empirical evidence on biological swarms reveal that interaction between individuals is selective. Each individual's neighbor is selected based on one or more featured factors. Based on the self-propelled model, we develop a general probability neighbor selection framework to study the effect of four typical featured factors (i.e., distance, bearing, orientation change and bearing change). In this work, two common cases are involved to comprehensively analyze the impact of the four featured factors on the collective motion. One is the flocking, the other is the responsivity to stimulus. The impact of different selection strengths of the featured factors on both cases are investigated. The effect of noise on flocking and different stimulus intensities on responsivity to stimulus are analyzed. This study allows us to get the insight of selective interaction and suggests the potential solution to overcome the trade-off between flocking and responsivity quality.
Collective motion is a common phenomenon in nature, such as the colony of bacterial clusters 1 , school of fish 2 , flock of birds 3 and active granular media 4 . Many efforts have been devoted to uncovering the mechanism of collective motion. Reynolds first proposed an approach to reproducing the flock of birds based on three simple rules (i.e., separation, alignment and cohesion), namely the Boid model 5 . Vicsek et al. 6 proposed the simplest self-propelled model to investigate the phase transition in the non-equilibrium system and reproduced the ordered movement by velocity alignment. Aldana et al. 7 described the original Vicsek model as the network system and found that the way that noise introduced into the system determines the order of phase transition in the self-propelled model. Chaté et al. 8 studied the Vicsek model with cohesion and formed the liquid-gas phase transitions. Couzin et al. 9 assumed the interaction between individuals is divided into three non-overlapping behavioral zones (i.e., zone of repulsion, orientation and attraction) and produced several complex collective behaviors by adjusting the width of each zone.
In recent years, plenty of empirical evidence has revealed selective interaction phenomena in animal groups [10][11][12][13] . Instead of considering all sensible neighbors, an individual selects its neighbors to react according to one or more featured factors. Such featured factors could be transient or dynamic determined by the relative motion of the focal individual. The transient featured factors are the stationary states of neighbors at a moment, such as the distance 5,6,9,14 and bearing [15][16][17] . On the other hand, the dynamic featured factors are usually coupled with the time, which suggests that neighbor selection depends on the change of the neighbors' states within a period, e.g., orientation change [18][19][20] and bearing change 21,22 . Although different featured factors have been observed, there is no agreement on how to use these factors when modeling a desired collective motion. Comprehensive analysis is necessary to obtain an insight view about selective interaction.
To this end, in this paper we firstly developed a general probability framework of neighbor selection to study the effect of four typical featured factors (i.e., distance, bearing, orientation change and bearing change) on the collective motion based on the self-propelled model. Then, we considered two common cases of collective motion, i.e., the flocking and responsivity to stimulus when comprehensively analyzing the effects of featured factors. Specifically, the flocking suggests that individuals form the highly parallel-group 9 from a disordered state. For the responsivity to stimulus, the stimulus can be interpreted as a kind of environmental cue 23 , which triggers the collective turn of the group. Additionally, the effects of different selection strength (i.e., the probability of being selected) of featured factors were analyzed in both cases. Also, different intensities of noise and stimulus were introduced when considering flocking and responsivity to stimulus, respectively. Interestingly, we found that solely increasing the strength of featured factors is harmful to the flocking; meanwhile, dynamic featured factors (i.e., orientation change and bearing change) are more responsive to the stimulus.
The rest of this paper is organized as follows. In "Methods" section, the general probability framework of neighbor selection is introduced to model self-propelled particle swarms. In "Results" section, the effects of each www.nature.com/scientificreports/ featured factor are analyzed from the perspective of selective interaction. The conclusions and further discussion are summarized in "Discussion" section.

Methods
Here we consider a group composed of N individuals moving in the continuous two-dimensional X-Y plane without boundary. The position of the i-th individual changes according to where v i is the velocity of the i-th individual. Each individual is self-propelled with a preferred speed v 0 and interacts with neighbors under random noise η (uniform distribution). The individual's sensing radius is an infinite range, which suggested that all individuals could be the neighbor of individual i. The velocity of the i-th individual can be calculated as follows: where γ 0 is the relaxation rate, which characterizes the self-propelled force. An individual will reach the preferred velocity v 0 faster with a larger γ 0 and vice versa. The second term F i accounts for the interactions between individuals. For simplicity, we assume that F i only depends on velocity alignment. Thus, F i is defined as K represents the strength of alignment and v d is the velocity with which v i should align.
Characterizing the neighbor selection according to different featured factors in a general way is the main part of our modeling framework. Inspired by Bode et al. 24 , we develop a probability framework. For a given individual j, the probability of being selected by the i-th individual obeys the following five selection models. Figure 1 shows the geometry of each featured factor.

Random selection model (RAND-S model).
In the RAND-S model, the focal individual selects one neighbor randomly (uniform distribution) at each time step. This model is actually equal to the Vicsek model (all sensible neighbors' states are averaged) if we consider time accumulation. In this paper, we use the RAND-S model as the baseline.

Distance-based selection model (DIST-S model).
In the DIST-S model, the focal individual selects one neighbor at each time step according to the relative distances between the focal and its neighbors. The distance-based selection rule was first proposed by Reynolds 5 and has been widely used to investigate the interaction network and the emergence of collective motion 6,9,14,25 . The probability of being selected in the DIST-S model is given as follows: where x i and x j are the position of individual i and j, respectively. c d ≥ 0 is a parameter of tuning selection strength on the nearest neighbor. When the c d is small, individuals prefer the distant neighbors. On the contrary, with a large c d , individuals tend to choose adjacent neighbors.

Bearing-based selection model (BEAR-S model).
In the BEAR-S model, the focal individual selects on neighbor at each time step according to the bearing between the focal and its neighbors. The BEAR-S model is inspired by the evidence of anisotropic sensory behavior in some animal species, which suggests that neighbors' bearing has an influence on the neighbor selection. Especially, the anisotropy of frontal preference has been (1) www.nature.com/scientificreports/ found to have significant effect on the collective motion. For example, Pita et al. 15 found that both the golden shiner and zebrafish have an acute vision in the fronto-dorsal region, which is in favor of group coherence. Lukeman et al. 17 found that the interaction of surf scoters is within the frontal sector of ± 30 • , which facilitates the collision avoidance in the group. The corresponding probability p bear ij is described as follows: where v i represents the velocity (unit vector) of individual i, and x ij is the unit position vector directed from i to j. ∠(·, ·) represents the angle between two unit vectors. The strength of attention in front can be adjusted by c b ≥ 0 . For a large strength c b , the focal individual tends to choose the neighbor that appeared straight ahead.
Dynamic-orientation selection model (DYO-S model). Inspired by the result 18 that the selection on neighbors is related to the orientation change of individuals rather than the distance in shoals of rummy-nose tetra species (Hemigrammus rhodostomus). Thus, we consider the orientation change within a period as a featured factor. For simplicity, we consider the time period as a single preceding step. The probability of DYO-S model is expressed as follows: where t is the step time interval in simulations. The v j t and v j t−1 are the velocity of individual j at time t and t − 1 . c o ≥ 0 is used to adjust the strength of selection. For a large c o , the focal individual tends to select the neighbor with the large orientation change.

Dynamic-bearing selection model (DYB-S model).
Much empirical research has revealed that the animal could be sensitive to the bearing change of their neighbors through the visual pathway 21,22 . The change in bearing is the manifestation of the apparent movement, which conveys information that can directly influence the individuals' movement. To be general, we assume that the time period of the bearing change is the same as that in the DYO-S model (i.e., one preceding step). The probability p dyb ij of selecting individual j can be described as follows: where x ij t and x ij t−1 are the position vector directed from i to j at time step t and t − 1 , respectively. c m is used to adjust the selection strength from small to large. When the c m is large, the focal individual tends to choose the neighbor with large bearing change.
To emphasize on the effect of the featured factors, we assumed that individuals only select one neighbor at each time step. Particularly, our framework includes the tunable parameter that controls the selection strength of the featured factors, which is able to select that specific neighbor from randomness to exclusiveness. In brief, our selection models comprise three steps.
(1) Calculation on the cumulative distribution function (CDF) of being selected based on a certain featured factor (Eqs. (3)-(6)), (2) Choose one neighbor by sampling from the calculated CDF, (3) Update the position and velocity of each individual using Eqs. (1) and (2).

Results
To investigate the property of each selection model, two cases are taken into account: flocking and responsivity to stimulus.
Flocking. It is common for animal groups to take off at the same time and form the ordered movement spontaneously, which is crucial to maximizing the chance of survival 26,27 . In the flocking case, we measure the ability to form the ordered movement from the random initial state. Here, we consider two experimental parameters: the intensity of noise η and the selection strength. The selection strength of each model is termed as follows: Two metrics are used to evaluate the simulation results. First, the polarization φ 6 is used to evaluate the degree of consensus, which represents the flocking quality. The polarization φ is defined as follows: here φ ∈ [0, 1] . When the group moves along with same direction, φ = 1 . Otherwise, if individuals diffuse with different directions, φ ≈ 0 . Second, we evaluate the flocking efficiency using convergence time T op 28 which is the time consumption before the polarization φ first reaches the given threshold. The definition of T op is given as follows: www.nature.com/scientificreports/ ε flock is the threshold of the ordered movement ( ε flock = 0.95 in this paper).

Responsivity to stimulus.
For biological swarm, the group constantly changes the moving direction to avoid the attack from predators. The responsivity to such sudden turning represents the effectiveness of information transfer within group 29 . To investigate the responsivity of featured factors, we consider the directional information as the stimulus. As suggested in the previous study 30 , a predefined number (i.e., IF num ) of informed individuals are randomly chosen from the group and abruptly turn with pi/2 respect to the current average moving direction. Moreover, we introduce stimulus intensity (i.e., the number of informed individual IF num ) and selection strength to conduct insightful analysis. Here, we use response accuracy δ group and response time T turn to evaluate the quality and efficiency of the responsivity to stimulus, respectively. The response accuracy δ group is defined as follows: is the average difference between the stimulus direction h s with the velocity of each individual v i . h s is orthogonal to the group's moving direction. A no represents the degree of consensus between the group's moving direction (averaged direction of all individuals) and the stimulus direction h s at one preceding step of the informed individual's turning moment, which is defined as . v no is the group's moving direction before the abruptly turning of the informed individual.
The time cost T turn during the response to stimulus is given as follows: here ε turn is the threshold if the group completes the response. t group is the turning moment of the informed individual. We set ε turn = 0.9 and t group = 10 in the following simulations.
Simulations are run with the following parameters: N = 100 , t = 0.02 , v 0 = 5 , γ 0 = 5 and K = 3 . The number of informed individual IF num is no more than 5% of total number ( IF num ∈ [1,5] ). All simulations are executed for 800 steps. Unless otherwise specified, any selection model discussed in this paper is the same as the RAND-S model if the selection strength is equal to zero.
Distance-based selection model. As for the DIST-S model, the preference on the nearest neighbor weakens the ability to form the ordered movement. Under the non-noise condition, when the selection strength is small ( c d < 1.5 ), the polarization φ and convergence time T op keep around the optimal value ( Fig. 2A-(1-2)), which suggests that individuals achieve the flocking efficiently ( Fig. 2A-(4)). However, when the selection strength is large ( c d > 1.5 ), the group splits into a few subgroups with different sizes and moving directions ( Fig. 2A-(3-5)), which implies the declination of flocking quality. Moreover, as the further increase of selection strength c d > 3 , individuals are scattered around and completely collapse ( Fig. 2A-(6)). This is because the focal individual overlooks the other individuals and only aligns with the nearest neighbor, which leads to the pairwise interaction. As a result of that, the slight directional deviation between individuals is hard to eliminate in time, resulting in the transform of group formation from coherence to disperse and the decay of flocking quality.
Under the noise condition, it is inevitable that the polarization and flocking efficiency decline with the noise increases. As shown in Fig. 2B-(1-2), the movement of individuals is disrupted by large noise intensity ( η > 2 ). When the noise intensity and selection strength is both small ( η < 1 and c d < 1.5 ), the group is able to form the ordered movement efficiently. However, when the selection strength c d > 1.5 , the polarization and flocking efficiency significantly reduce under the same noise intensity ( η < 1 ). To sum up, the strong preference on the nearest neighbor has a negative impact on noise tolerance.
The responsivity to stimulus is promoted by the preference on the distant neighbor rather than the nearest neighbor. According to Fig. 2C-(1), when the stimulus intensity is weak ( IF num = 1 ), the response accuracy δ group never exceeds 0.5, which means that the group is unable to response to the stimulus (Fig. 2C-(3)). Additionally, from Fig. 2C-(1-2), it can be seen that the response accuracy and efficiency show an apparent upward trend with the stimulus intensity increases, which implies that individuals require more stimulus information to make response. The trajectory of the response under a large stimulus intensity ( IF num = 5 ) is shown in Fig. 2C-(4). However, under the maximum stimulus intensity ( IF num = 5 ), the response accuracy and efficiency still decrease consistently as the increase of c d . Basically, the reason for low responsivity is that individuals only focus on their nearest neighbors during the response process, which blocks the stimulus information transfer within the group. In other words, directional information is hard to efficiently spread in the DIST-S model.
Bearing-based selection model. For the BEAR-S model, the frontal preference has a negative impact on forming ordered movement. In the non-noise situation, according to Fig. 4A-(1-2), when the frontal preference is weak ( c b < 10 ), individuals achieve the ordered movement efficiently. As the increase of c b , the group tends to split into two or more clusters. These clusters are still likely to form the ordered movement when the frontal preference is at an intermediate level ( c b ≈ 20 ), as shown in Fig. 4B-(1). When the selection strength is large ( c b > 25 ), the group splits into two clusters with totally different moving directions (Fig. 4B-(2)). Meanwhile, the polarization φ declines to 0.6, and the convergence time T op rises to 600 simulation steps. It arises from the fact that the frontal preference restrains the individuals' interaction range. The neighbor in front directly affects www.nature.com/scientificreports/ the individual at the back. On the contrary, the back-end individual never has an influence on the neighbor in front. In other words, the interaction becomes regional and its topology switches slowly. This leads to the group dispersion and it is difficult for individuals to keep ordered. In the noise situation, it is interesting that different strength of frontal preference leads to quite different decay tendencies. When the selection strength is small ( c b < 10 ), the polarization and flocking efficiency decline monotonically with the increase of noise intensity. Counterintuitively, under the strong frontal preference ( c b > 10 ), the BEAR-S model exhibits a non-monotonic decline tendency. Specifically, as shown in Fig. 4A-(3-4), when the noise intensity varies in η ∈ [0, 1] , the polarization φ takes an upward trend while the convergence time T op has a rapid decrease. This suggests that under this situation the flocking quality and efficiency have been promoted evidently. The polarization and flocking efficiency achieve the optimal value at η ≈ 0.8 rather than that under the non-noise condition ( η = 0 ). It is because that the involvement of noise makes it easier for individuals to interact with more neighbors. As shown in Fig. 4B-(3-4), the neighbor density under the noise η = 0.1 (Fig. 4B-(3)) is narrower than that under the large noise intensity η = 1 (Fig. 4B-(4)). The individuals' interaction range is expanded to the larger level, which means that the flocking quality is promoted as the result of noise perturbation.
Frontal preference is negative for the stimulus transfer within the group. According to the Fig. 4C-(1-2), with the increase of c b , the general trend in response accuracy δ group decreases consistently while the response time T turn shows an apparent growth for any number of the informed individual. More evidence can be seen in Fig. 4C-(3), individuals with strong frontal preference barely respond to the stimulus. The group cannot follow the turning of the informed individuals until the stimulus intensity increasing to IF num = 5 (Fig. 4C-(4)). Nevertheless, with the increase of c b , the response accuracy and efficiency still decrease rapidly. It is due to that the strong frontal preference limits the individuals' perception and breaks the stimulus propagation link within the A B C (1) (2) (1) (3) ( 4) www.nature.com/scientificreports/  www.nature.com/scientificreports/ group. As shown in the Fig. 4C-(1-2), the variance of response accuracy δ group and response time T turn are both high. There is another interesting finding shown in Fig. 3, the informed individual in the front of the group does promote the responsivity to stimulus, but it is still unable to transfer the stimulus effectively on occasions. The reason might be that the frontal preference is highly related to the neighbors' bearing. Thus, the group responsivity becomes sensitive to the the initial position of stimulus (i.e., the spatial distribution of the informed individual).  ≈ 60 ), the polarization φ declines to 0.65. At the same time, the group is slightly dispersed and keeps the ordered formation ( Fig. 5B-(1)). As for the flocking efficiency, as shown in Fig. 5A-(2), the convergence time T op rises consistently with the increase of c o . The decrease of flocking efficiency is due to that the orientation change of each individual is lack of discrepancy. The switching of neighbors is driven by the stochastic factor in simulations. The effectiveness of neighbor selection is heavily reduced, which leads to the blind movement.

Dynamic-orientation selection model.
In the noise situation, the noise tolerance is impaired by the large selection strength c o . According to Fig. 5A-(3), when the selection strength c o is small ( c o < 1 ), the group is able to keep the high polarization even with the Neighbor Density Neighbor Density www.nature.com/scientificreports/ noise intensity η = 1.5 . However, with the increase of c o , individuals become vulnerable to the noise. As for the polarization, the decline tendency shifts from a gradual decrease to a rapid drop as the approximately discontinuous transition manner. The critical noise intensity of the disorder-order transition is decayed from η ≈ 1.5 ( c o = 1) to η ≈ 0.8 ( c o = 40). The DYO-S model's intrinsic mechanism (i.e., identification of the neighbor with a large orientation change) is responsible for the significant impairment of the noise tolerance. Individuals suffer from the false information (i.e., the wrong alignment direction) under the large selection strength c o , which means that the group cannot form the ordered movement ( Fig. 5B-(2)). As for the flocking efficiency, as shown in Fig. 5A-(4), with the increase of noise, different selection strength c o results in the different decline tendencies of flocking efficiency. Specifically, when the selection strength is small ( c o < 1 ), it is reasonable to find that the convergence time T op maintains an upward trend with the increase of noise intensity. Counterintuitively, with the large selection strength ( c o ≈ 40 ), the effect of noise on the convergence time T op shows a non-monotonic decline tendency. Interestingly, as shown in Fig. 5A-(4), there is an unexpected decrease in the convergence time T op with the noise intensity ( η ∈ [0, 0.5] ). Such promotion of flocking efficiency is mainly attributed to that the neighbor switch becomes more frequent under the noise condition, which implicitly expands the interaction range. As shown in Fig. 5B-(3-4), comparing with the condition of noise intensity η = 0.1 , the interaction range is expanded by larger noise intensity η = 0.5.
The large selection strength c o improves the response accuracy and efficiency to the stimulus. As shown in Fig. 5C-(1-2), under the low stimulus intensity ( IF num = 1 ) and also the increase of c o , the response accuracy δ group rises up fast while the response time T turn decreases significantly. Even if the selection strength is large ( c o > 10 ), the response accuracy δ group and time cost T turn remain stable at their optimal values regardless of stimulus intensity changes. However, when the selection strength c o < 10 , the response accuracy and efficiency have two different trends depending on the different levels of stimulus intensity. On the one hand, when the stimulus intensity IF num ≤ 2 , the response accuracy and efficiency increase monotonically to the optimal value. On the other hand, with the large stimulus intensity ( IF num ≥ 3 ), the promotion of responsivity shows a nonmonotonic tendency (i.e., first decreasing and then increasing). For example, when the stimulus intensity IF num is up to 3, it is sufficient for the group to respond under the influence of average interaction ( c o = 0 ), as shown in Fig. 5C-(3). Since the focal individual takes the limited focus on the neighbor with large orientation change ( c o < 5 ), the influence of average interaction is neutralized and leads to the reduction on responsivity to stimulus ( Fig. 5C-(4)).
Dynamic-bearing selection model. For the DYB-S model, the large selection strength c m leads to the decrease in polarization and has an unfavorable effect on the flocking efficiency. In the non-noise situation, as shown in Fig. 6A-(1), with the increase of c m , the polarization falls gradually at φ ≈ 0.8 . Meanwhile, as shown in Fig. 6A-(3-7), the group disperses laterally as a line-like formation but still forms the ordered movement. According to Fig. 6A-(2), the convergence time T op shows an evident upward trend. When the selection strength c m > 20 , the flocking efficiency reaches the lowest level ( T turn ≈ 500 simulation steps). The reason for the decrease in flocking efficiency is that preference on neighbors is changed. Specifically, when the selection strength c m is large, the individuals' preference on neighbors converts from the bearing change to the close proximity. Additionally, when traveling with the same distance, the bearing change is more noticeable for the closer individuals, as shown in Fig. 6A-(8). This conversion leads to the regional interaction between individuals. It is difficult for the group to keep coherent and achieve the ordered movement efficiently.
In the noise situation, the group that pays more attention to the neighbor with the large bearing change has less tolerance on the noise. From Fig. 6B-(1), we can see that the polarization declines from gradually to rapidly with the increase of c m . When the selection strength c m is small ( c m < 2 ), individuals are able to keep the high polarization and flocking efficiency when the noise intensity η < 1 (Fig. 6B-(1-2)). When the selection strength c m is large ( c m > 2 ), polarization and flocking efficiency declines accordingly under the same noise intensity ( η < 1 ). This is because that the noise affects the neighbors' bearing perception and confuses individuals during the neighbor selection process.
The preference on the neighbor with the large bearing change allows the group to transfer the stimulus accurately and efficiently. According to the Fig. 6C-(1-2), even with the low stimulus intensity ( IF num = 1 ), the response accuracy δ group increases rapidly and the response time T turn falls sharply as the increase of c m . When the selection strength c m > 2 , the response accuracy δ group and efficiency reaches their optimal value, respectively. Moreover, with the increase of stimulus intensity, the responsivity of the group keeps at the optimal level regardless of any selection strength c m ( c m > 0 ). The reason for such great responsivity is that the interaction topology that emerged in the group facilitates the stimulus transfer. In other words, the interaction topology based on bearing change can efficiently spread the directional information within the group.

Discussion
In summary, we analyzed the effect of different featured factors on the collective motion based on the probability framework, as it is necessary to get the insight of the selective interaction.
For the flocking, in the non-noise situation, it is important to find that the increase of selection strength has a negative impact on the flocking for both transient and dynamic featured factors. In the noise situation, it is interesting that the proper noise intensity promotes the flocking quality in both the BEAR-S model and the DYO-S model.
On the other hand, the responsivity of different featured factors differs greatly. The information transfer by dynamic featured factors (i.e., orientation change, bearing change) is more efficient than that through the transient. Generally, the most significant difference between the dynamic and transient featured factors is the introduction of the change of neighbor's state. The dynamic featured factors consider the neighbors' state at not www.nature.com/scientificreports/ only a particular moment but also the collection of the neighbor's state over a period. The dynamic featured factor represents the gradient of the stimulus. Comparing with the transient featured factors, the dynamic featured factors are intrinsically suitable to transfer the stimulus within the group. According to the simulation results, we found that there is a trade-off between the flocking quality φ (i.e., persistence) and the responsivity quality δ group . Especially for the DYO-S and DYB-S model, the group is bound to choose either high flocking quality or high responsivity when they exhibit collective motion. Persistence and responsivity are two distinct perspectives of collective motion, which can not be achieved simultaneously. The persistent (coordinated) movement is achieved by ruling out deviations from group's moving direction, which demands the overdamped interaction dynamics. However, the responsivity is obtained by spreading fast turns across the group, which requires the underdamped interaction dynamics 23,31 . In this work, our results suggest that there may be two possible ways to reconcile this trade-off in the view of selective interaction. For one thing, the combination of different featured factors could be an effective solution. As an instance, the DYB-S model (1) (1) (3) (4) www.nature.com/scientificreports/ incorporated with the proper frontal preference (i.e., BEAR-S model) is likely to promote the noise tolerance and remain the high responsivity to stimulus. For another thing, featured factors shift in accordance with different tasks. At the beginning of the flocking, it is reasonable to adopt the distance in the collective model for its simplicity. In addition, the neighbor's bearing (i.e., the frontal preference) could be involved to improve the noise tolerance in the complex environment. Importantly, when individuals are under stimulus, the featured factors of change in orientation and bearing have potential to promote the responsivity of the group. Switching the featured factors in the collective model will make it easy to handle different situations.
In this paper, we only involve the single-neighbor interaction in our proposed framework. However, the number of selected neighbors is another significant part of the neighbor selection strategy. It has been shown that the emergent coordinated movement depends strongly on how many neighbors a focal individual can pay attention to in Refs. [32][33][34] . To be general, we also conduct some preliminary experiments on the interaction with multiple neighbors in the case of flocking. The main conclusions remain consistent with the experiments of single-neighbor interaction. The general trend of flocking quality φ and efficiency T op declines with the increase of selection strength for any number of selected neighbors. Interestingly, with the increase of the number of selected neighbors, the declination trend of flocking quality and efficiency becomes much smoother, which suggests that further analysis is necessary to have a better understanding on the effects of multi-neighbor interaction. Additionally, there is no agreement on the number of selected neighbors in the biological swarm, which might be different from animal species, e.g., the starling birds usually interact with 6-7 neighbors 16 , while the interaction in the fish group appears to be restricted to a rather low number of neighbors (one or two individuals) 14,18 . Therefore, we plan to comprehensively investigate the influence of the number of selected neighbors on a robust ordered motion in the future.

Data availability
The datasets generated and/or analysed during the current study are available from the corresponding author on reasonable request.